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Abstract 

A ID metallic surface state was created on an anisotropic InSb(OOl) surface covered with Bi. 
Angle-resolved photoelectron spectroscopy (ARPES) showed a ID Fermi contour with almost no 2D 
distortion. Close to the Fermi level (Fp), the angle-integrated photoelectron spectra showed power- 
law scaling with the binding energy and temperature. The ARPES plot above E-p obtained thanks 
to thermally broadened Fermi edge at room temperature showed a ID state with continuous metallic 
dispersion across Ep and power-law intensity suppression around Ep. These results strongly suggest 
a Tomonaga-Luttinger liquid on the Bi/InSb(001) surface. 

PACS numbers: 71.20.-b, 73.20.At, 79.60.-i, 71.10.Pm 
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The Fermi liquid theory of ordinal three-dimensional (3D) metals breaks down in one¬ 
dimensional (ID) systems to produce various exotic quantum phases. Tomonaga-Luttinger 
liquid (TLL) [l, 2| is an exactly solvable model of a gapless ID quantum system, that 
is characterized by power-law scaling and spin-charge separation for low-energy excitation 
spectra j^. So far, various ID systems have been studied as TLL candidates. Only a few 
have shown metallic states with power-law spectral features, e.g., carbon nanotubes ^ and 


lithium purple bronze 


js, 6|. The other characteristic of spin-charge separation has been 


reported in some 3D materials with ID electronic structures [7 h9(]. 

The surface of semiconductors is known to show various self-assembled ID atomic struc¬ 


tures that are reg arded as suitab 


e systems for studying ID metallic states 


QQ, 


such 


as In/Si(lll) 12|, Au/Si(557) 13|, and Pt/Ge(001) 14|. On such surface systems, in-situ 


electron/hole doping can be performed by the deposition of additional atoms, and the local 
atomic structure can be observed/controlled by scanning probe techniques. These manipu¬ 
lations on the surface ID states can provide further insight into ID physics. Hence, surface 
TLLs have attracted much attention in the past two decades. However, despite the various 
ID surfaces, most of them do not behave as TLLs at low temperature. Possible reasons for 
the non-TLL behavior have been discussed for each specific system, such as a metal-insulator 
transition at low temperatures and interaction with other two-dimensional (2D) surface 


bands 


16|. 


So far, the only candidate for a surface TLL is an Au nanowire assembled on Ge(OOl) 
(Au/Ge(001)), which shows power-law scaling for low-energy spectral features according 
to scanning tunneling spectroscopy (STS) and angle-resolved photoelectron spectroscopy 
(ARPES) [l7|, ll8|. However, a recent work presented doubt with regard to the TLL picture 
and proposed an alternative interpretation of the STS spectra: two different and gapped 
states lying just above/below the Fermi level {Ep) [l9|. This picture can also explain the STS 
spectral shape around Ep. In order to distinguish these two possibilities without ambiguity, 
the ID band dispersion needs to be traced across Ep continuously with k resolution in 
reciprocal space. However, to the best of our knowledge, no such experimental data have 
been reported about surface TLLs so far. While some ARPES results have been given for 

nnn 

Au/Ge(001) [17|, Il8|, 12^ , the dispersion above Ep is not accessible with these data. 

In this letter, we report on the case of a new surface TLL candidate discovered 
with ARPES on a Bi-induced anisotropic structure on InSb(OOl) (Bi/InSb(001)). The 


2 






























Bi/InSb(001) surface state showed a ID Fermi contour (FC) with almost no undulation 
and stayed metallic down to 35 K. The ID feature is the unique contour on Ey of the 
sample, which showed a parabolic dispersion with its bottom at the center of the surface 
Brillouin zone (SBZ). The photoelectron spectra in the proximity of obey power-law scal¬ 
ing as a function of binding energy and sample temperature with the same power index a = 
0.7 ±0.1. The ARPES plot divided by the Fermi distribution function at room temperature 
(RT) indicated the ID state with continuous metallic dispersion across Ey and power-law in¬ 
tensity suppression around Ey- These results strongly suggest that the ID surface structure 
of Bi/InSb(001) hosts a surface TLL. 

The Bi/InSb(001) surface was prepared by evaporating three monolayers of Bi on the 
clean surface of InSb(OOl) substrates (n-type, carrier concentration of 2x10^^ cm“^) and a 
subsequent flash by direct current heating up to 680 K for ~10 seconds. After the prepa¬ 
ration, the p(lx3) surface periodicity was observed by low-energy electron diffraction (Fig. 
1(a)). The detailed procedure on sample preparation and characterization, especially the 


comparison with the c(2x6) surface reported in a previous work [l|, is shown in the supple¬ 
mental material (SM) 22|. 

Figs. 1(c)(e) indicate the surface electronic structure of Bi/InSb(001) measured with 
ARPES. ARPES measurements were performed with a He lamp and synchrotron radiations 
at the CASSIOPEE beamline of synchrotron SOLEIL and BL7U of UVSOR-III. The photon 
energies used ranged from 12 to 40.8 eV. The photoelectron kinetic energy at Ey and the 
overall energy resolution of each ARPES setup were carefully calibrated by the Fermi edge of 
the photoelectron spectra from freshly evaporated gold on the sample after the measurement 
(see Fig. 3(a)). 

Figure 1 (c) shows the Fermi contour (FC) around Ey as measured with linearly polarized 
photons, kx and ky are dehned to be parallel to [110] and [110], respectively. The linear 
shape of the FC along ky clearly indicated a ID metallic state on the surface. Such a ID FC is 
consistent with the anisotropic (1x3) surface periodicity indicated by FEED. The dispersion 
of the ID state did not change with the photon energy (shown in SM 22|), which indicates 
no dispersion along k^ as would be expected for a surface state. Figure 1 (c) also shows that 
there were no other metallic states on the surface. Therefore, this FC provides a unique 
and highly anisotropic conduction path on the surface. As shown in Fig. 1(d), the ARPES 
momentum distribution curves (MDCs) indicated that the SS was exactly ID without any 
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2D distortion. MDCs were captured in the second SBZ because the SS dispersion became 
clearer there owing to the photoexitation matrix element effect. The lack of 2D distortion 
can also be observed in the E-ky plot in Fig. 1(e). The dotted curve is the upper edge of the 


projected bulk valence bands as calculated according to Ref. 


24l |. The bulk valence band 


maximum position was determined with the linear method (see SM for details 22|). 

Note that the ID FC ap peared with a single domain, in contrast to the double-domain 


case of Au/Ge(001) 18|, [2^, where the direction of the conduction path is stil 


because of the two equivalent domains rotated 90° with each other 


19 


25 


under debate 


261. This is 


because of the reduced symmetry of the zinc-blend crystal structure compared to group-IV 
semiconductors. Along the [001] direction, a III-V semiconductor lattice is stacked alter¬ 
nately with group III and V atoms (indium and antimony for InSb). Hence, the equivalent 
plane appears at every two atomic steps on the (001) surface of such crystal. In the group 
IV semiconductor crystal, every atomic plane along [001] is equivalent and appears with 90° 
rotation at each atomic step to form the double-domain surface structure. The 90° rotation 
at each atomic step also occurs for III-V semiconductors; the double atomic step results in 
a 180° rotation, so the steps are actually identical thanks to the two-fold rotation symme¬ 
try. Thus, the conduction path of the Bi/InSb(001) surface was clearly demonstrated to be 
parallel to [110]. 

Figure 2 shows the dispersion of the surface band on Bi/InSb(001). The surface state 
formed a single parabolic band with its bottom at fca, = 0 A“^. This band is the unique 
metallic feature on the surface. This ID metallic band crossed E-p only once on each side of 
the SBZ. At higher binding energies, the bottom of the surface-state band showed a small 
shift (~50 meV) between ky = 0.0 and 0.2 A“^. This would be due to an overlap of the 
bulk valence bands. The dotted (solid) curves in Fig. 2 are the upper (lower) edges of the 
projected bulk valence (conduction) bands and were calculated in the same manner as those 
in Fig. 1(e). The bottom of the SS overlapped the bulk valence bands at the center of the 
SBZ. Note that the bulk conduction band does not cross Ep, as shown by the solid curve in 
Fig. 2(a). 

In order to examine the ID nature of the surface state, we compared the angle-integrated 
photoelectron spectrum with that from ordinary metal (An evaporated on the sample after 
the ARPES measurement), as shown in Fig. 3 (a). The angle integration was performed 
at ky = 0.25±0.05 A“^. Figure 3 (a) clearly shows the suppression of the photoelectron 
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intensity at E-p compared to that from An, which suggests a deviation of the surface state 
from the Fermi-liquid framework. We tried to fit the spectrum to the TLL spectral function 
considering the hnite temperature, as given in some of the literature j^, 28|. Based on ref. 


27|, the spectral function of the angle-integrated PES can be written as, 


I{e,T) (xT^cosh (^0 


1 + a . e 
2 


/(e,T), 


( 1 ) 


where e = E/k-pT is the temperature-normalized energy, P is the gamma function, and 
f{e,T) = (e^ -|- 1)“^ is the Fermi distribution function. When convolved with the Gaussian 
instrumental energy resolution {AE = 10 meV), this function reproduced the spectra well 
at 35 K (Fig. 3 (a)) with a = 0.65 ± 0.05 in the low-energy region. This function only 
hts the spectrum below 0.1 eV. This was because the TLL framework is only valid in the 
proximity of Ep (see SM for detailed analysis |22| ). 

The power-law scaling of a TLL should show an universal power index with various 
parameters. Figs. 3 (b) and (c) show the temperature dependence of the photoelectron 
intensities at Ep. The intensity from 180 to 50 K decreased linearly with log T, which 
indicates power-law scaling of the spectral intensities, Ie=Ef(T) oc T". The slight change in 
the slope below 50 K was due to the finite AE. This broadening was already convolved in 
the htting curve (blue solid line) to give a power index of a = 0.71 ± 0.10, which is a value 
very close to that obtained by fitting the energy scale (Fig. 3 (a)). Note that Ie.eAT) 
should be constant in the case of a Fermi-Dirac system. 

Another method to examine universal power-law scaling is to normalize the spectra taken 
at various temperatures to T“ and plot them versus the temperature-normalized energy e, as 
shown in Fig. 3 (d). With such a renormalization, all spectra with various T become almost 
identical in the low-energy region, typically in e < 5. For such identification, the scaling 
actor a = 0.75 ± 0.05 is required (spectra renormalized with another a are shown in the SM 


22l|L Because such a temperature normalization method can include the intensity not only 


at E=Ep but within the finite energy window, this renormalized plot would be the best way 
to estimate the power index from angle-integrated spectra taken at various temperatures. 
All of these spectral features indicate that the ID surface state on Bi/InSb(001) obeys an 
universal power-law scaling with a = 0.7±0.1, which agrees with that expected for a TLL. 

To obtain further insight into the ID SS, we looked into the detailed structure of ARPES 
spectra. Figure 4 (a) shows the ARPES intensity plot at RT along kx divided by the Fermi 
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distribution function convolved with the instrumental resolution. For better statistics, the 
ARPES plot was integrated along ky within \ky\ <0.3 A“^. The momentum distribution 
curves (MDCs) are shown in Fig. 4 (b). They are htted well by a pair of Voigt peaks 
(Gaussian and Lorentzian widths are hxed for every peak) and linear backgrounds (solid 
lines in Fig. 4 (b)). Weak residuals around = 0 at the binding energies below -0.1 
eV would be due to the bulk conduction bands. From the MDCs, we obtained the peak 
positions at each binding energy and overlaid them on Fig. 4 (a). They indicated the Fermi 
wavevector |/cf| to be 0.05 ±0.02 A“^. This estimation is consistent with the ARPES data 
taken at low temperature shown in Fig. 2. This shows the continuous, metallic dispersion 
of the ID SS across Ep. Moreover, the ARPES intensity plot clearly shows that the spectral 
intensity of the ID SS was suppressed around Ep, without breaking the continuity of the 
band dispersion. Such suppression of the density of states (DOS) around Ep without any 
gap opening reflects the power-law scaling of the TLL. This result excludes alternative 
possibilities without ambiguity. Two different states below/above E-p cannot reproduce such 
a dispersion. In addition, the angle-integrated spectrum showed a temperature-broadened 
power-law shape for the density of states around Ep, as expected in theory 0 . At RT, 
it sometimes becomes difficult to determine the exact position of Ep because of the broad 
Fermi edge. However, note that such ambiguity did not affect the DOS suppression obtained 
by ARPES: the intensity suppression appeared irrespective of slight artihcial shifts in Ep. 

In summary, a new candidate for a surface TLL formed on Bi/InSb(001) has been dis¬ 
covered. A surface state formed a ID FC with almost no 2D undulation and stayed metallic 
down to 35 K. This ID FC is the only metallic feature of the sample and hence provides an 
unique ID conduction path on the surface. The photoelectron spectra in the proximity of Ep 
obeyed the power-law scaling as a function of binding energy and sample temperature with 
the same power index a = 0.7 ±0.1. Moreover, the ARPES plot above Ep obtained thanks 
to thermally broadened Fermi edge at RT indicated a ID state with continuous metallic 
dispersion across Ep and power-law intensity suppression around Ep. These results strongly 
suggest that the ID surface structure Bi/InSb(001) hosts a surface TLL. The new ID SS 
on Bi/InSb(001) would provide a fertile playground for further studies of low-dimensional 
physics. 
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Supplementary material for: Surface Tomonaga-Luttinger liquid state on 

Bi/InSb(001) 

Sample growth and characterization of Bi/InSb(001)-p(lx3) and c(2x6) surfaces 

InSb(OOl) substrates (n-type, carrier concentration of 2x10^^ cm“^) were cleaned by 
repeated cycles of sputtering and annealing up to 680 K until the c(8x2) low-energy electron 
diffraction (LEED) pattern was observed. Then, a nominal 3 monolayers (ML) of Bi were 
evaporated at room temperature: 1 ML is dehned as the atom density of bulk-truncated 
InSb(OOl) and the coverage was estimated by using a quartz microbalance. A subsequent 
flash by direct current heating up to 680 K for ~10 seconds produced the p(lx3) surface, as 
indicated by the sharp and low-background LEED pattern shown in Fig. 1 (a) of the main 
text. As shown in Fig. SI (a), the excess Bi was removed after the flash. From the height 
of Bi 5d with linear-background subtraction, the coverage of Bi on the (1x3) surface was 
~1.5 ML. 

Sometimes, we could hnd faint fractional spots showing the c(2x6) periodicity on the 
flashed surface, which suggested a small fraction of the coexisting area with the c(2x6) 
periodicity, as shown in Fig. SI (b). The previous work [l| reported a Bi/InSb(001)- 
c(2x6) surface prepared by 2 ML of Bi evaporation followed by annealing at ~ 470 K for 2 
hours. This surface was reported to be metallic based on the angle-integrated photoelectron 
spectrum taken at room temperature (RT). For comparison, we prepared the c(2x6) surface 
according to Ref. [l|. The LEED pattern showed the c(2x6) periodicity (see Fig. SI (c)). 
As shown in Fig. S2 (a), the Fermi contour measured by ARPES did not show any one¬ 
dimensional SS. Instead, there was a two-dimensional metallic state. The angle-integrated 
spectra in Fig. S2 (b) qualitatively agrees with that in the previous work. 

Based on these results, we concluded that the ID SS observed in this work comes from 
the flashed (1x3) surface and that the c(2x6) area does not contribute at all to the ID 
states. 

ARPES images taken with different photon energies 

Figure S3 shows the ARPES images taken with various photon energies. The parabolic 
state below the Fermi level showed no dispersion; thus, it is the surface state. 
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Bulk valence band maximum (VBM) 


Figure S4 shows the angle-integrated photoelectron spectra taken with four different 
photon energies. We determined the energy position of VBM to be the intersection between 
the linear hts of the background and the linear portion of the valence band leading edge 
[^. Based on the so-called linear method, the valence band maximum was estimated to be 
0.19±0.01 eV. 

Fitting cnrve for angle-integrated photoelectron spectra 

We checked the dependence of the power index to the energy window as shown in Fig. 
S5. Fitting was performed based on Eq. (1) in the main text. Below 0.1 eV, the power 
index was almost the same and the htting curve agreed well with the spectrum. With higher 
binding energies, the power index tended to become smaller and Eq. (1) failed to reproduce 
the experimental data. This must be because the TLL framework is only valid in the low- 
energy region. In addition, this deviation from the theoretical curve would be because the 
bottom of the SS was ~ 0.25 eV: around the bottom of the surface-state band, the dispersion 
is no longer linear. 

Renormalized photoelectron spectra with various scaling factors 

Figure S6 shows the universal scaling plot of the angle-integrated photoelectron spectra 
with different scaling factors a. Fitting of the spectra taken at various temperatures to the 
universal curve was optimized at a = 0.75. 
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FIG. 1: (Color online), (a) Low-energy electron diffraction patterns of Bi/InSb(001)-(lx3) taken 
at room temperature (RT). (b) Surface and bulk Brillouin zones of Bi/InSb(001). (c-e) ARPES 
results at RT measured with hv = 21 eV. (c) Constant energy contour around the Fermi level {E-p). 
The dashed lines indicate boundaries of the surface Brillouin zone (SBZ). kx and ky are defined 
parallel and perpendicular, respectively, to [110]. ARPES intensities are integrated over 20 meV 
windows centered at Ep. (d) Momentum distribution curves in second SBZ. (e) Band dispersion 
around Ep along ky. 


11 












































0.2 



FIG. 2: (Color online) (a, b) Band dispersions along at 8 K with hv = 15 eV. The dotted (solid) 
curves represent the upper (lower) edges of projected bulk valence (conduction) bands. 
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Nomialized energy (eW/k^J) 


FIG. 3: (Color online), (a) Photoelectron spectrum (red markers) integrated along kx {ky = 
0.25±0.05 A“^). The solid line between ± 0.1 eV is the fitting curve (see text for details), (b) 
Temperature dependence of the photoelectron intensity at Ep {Ie=Ef{T)). Each spectrum was 
normalized at the intensity at 0.2 eV below Ep in order to compensate the sample degradation. 
The solid (blue) line is the fitting curve based on power-law scaling (T^'^^) convolved with AE (10 
meV). The dashed line is without broadening from AE. (c) The same as (b) plotted in a log 
scale, (d) Universal scaling plot of the spectra at various temperatures, using a T-renormalized 
energy scale {e\/kpT). The thick dotted (green) line is from eq. (1) with the scaling factor a = 
0.75. 


13 




























(a) A, from SBZ center (A ') (b) from SBZ center (A') 



FIG. 4: (Color online), (a) ARPES intensity plot along kx divided by the Fermi distribution 
function convolved with the instrumental resolution taken at hu = 21 eV. ARPES data were 
integrated along ky within \ky\ < 0.3 for better statistics, (b) ARPES MDGs along kx within 
20 meV energy window. The intensities of each curve were normalized by the peak height. Solid 
lines are the fitting curves (see text for details), (c) An angle-integrated EDC within \kx\ < 0.3 
divided by the Fermi distribution function, similar to (a). 
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FIG. SI: (a) Core-level spectrum for the as-deposited and flashed Bi/InSb(001) surfaces taken at 
hu = 40.8 eV (He II). (b) FEED pattern of a flashed Bi/InSb(001) surface with coexisting c(2x6) 
spots (indicated by red markers). The pattern is distorted because of the flat MCP screen, (c) 
LEED pattern of the Bi/InSb(00I)-c(2x6) surface. 


-0.4 -0.2 0.0 0.2 0.4 



EIG. S2: ARPES results of the Bi/InSb(001)-c(2x6) surface with hu = 19 eV at room temperature, 
(a) Eermi contour, (b) Angle-integrated spectrum kx {ky = 0 A“^). 
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FIG. S3: ARPES intensity maps of Bi/InSb(001) along kx {ky = 0.2 A with various photon 
energies taken at 8 K. 



_l_I_L._I -- - ^ 

0.4 0.2 0.0 


Binding energy (eV) 

FIG. S4: Angle-integrated photoelectron spectrum along kx {ky = 0 A“^) measured at 13 K. The 
dashed lines are linear fits of the leading edges from the bulk valence bands. 
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FIG. S5: Angle-integrated photoelectron spectrum (same data as figure 3 (a)) fitted with various 
energy ranges. 
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FIG. S6: Universal scaling plot of the spectra at various temperatures with various scaling factors 
a, using a T-renormalized energy scale (eV/ZcsF). (a) a = 0.65, (b) a = 0.70, (c) a = 0.75, (d) a 
= 0.80, and (e) a = 0.85. 
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